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1. Introduction 

Loop Quantum Gravity (LQG) quantises canonical General Relativity, in a background 
independent and non-perturbative manner. It does so by considering the three- 
dimensional spatial slice to be a network of triads and then formulating the classical 
theory in terms of holonomies of the triad connections, around closed edges of the 
network, and fluxes of the triads, through the surfaces enclosed by these edges. Whilst 
the full, inhomogeneous theory has yet to be fully developed, significant progress has 
been made by applying the approach to symmetry reduced, mini-superspace models. 
In particular, by applying these principles of quantisation to Cosmology, one arrives at 
Loop Quantum Cosmology (LQC), which has proved itself especially successful |T|. 

Within LQC it has been possible, on the one hand, to explicitly show the 
avoidance of singularities that typically plague the classical versions of these cosmological 
theories |2l |3] and, on the other hand, to develop powerful effective theories, that capture 
first order corrections to classical cosmology In particular, it has been shown that, 
in order for the cosmological theory to be theoretically consistent |5l [6], the physical 
volume (rather than the scale factor, or the triad component) must be the fundamental 
variable used. This conclusion has been backed up by various phenomenological [3 IH] 
and theoretical [9] arguments. However, it is important to realise that since there is 
still lacking a complete derivation of LQC from the fundamental inhomogeneous and 
anisotropic theory, it is important to test LQC predictions for robustness (see e.g., 
Ref. do]). 

Various aspects of anisotropic cosmologies have been studied within LQC in the 
past [m [121 US]) however the first full and consistent quantisation of a Bianchi I 
cosmology (the simplest of anisotropic cosmological models) was achieved in Ref. [11] 0- 
Even more encouraging, the link back to the underlying full theory (LQG) has been 
strengthened, by considering the flux of the triads through surfaces consistent with the 
Bianchi I anisotropic case. With the quantisation of the Bianchi I model under control, 
it is possible to ask whether LQC features, obtained within the context of isotropic 
Friedmann-Lemaitre-Roberston- Walker (FLRW) cosmology, are robust, at least with 
respect to this limited extension of the symmetries of the system. 

In this paper, we will show that one can choose among different isotropic 
embeddings of the anisotropic Bianchi I model. We will then demonstrate that the choice 
of different embeddings has important consequences for the precise form of discretisation, 
produced by the loop quantisation procedure, in the isotropic sub-system. Note that, 
we use Bianchi I as an example of an anisotropic model. In other words, we do not 
imply that it is physically natural, nor that our results are constrained to hold only for 
this type of an anisotropic model. 

In Section 2, the basics of LQC will be sketched, for both isotropic and anisotropic 

I Performing a von Neumann stability analysis of the Hamiltonian constraint equation valid for 
anisotropic Bianchi I model, we have shown |15] that the the difference equation is unconditionally 
unstable. 
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models. In Section [31 it will be demonstrated that an isotropic embedding of an 
anisotropic system can be realised in various ways, leading to consequences for the 
form of the induced isotropic system. This will be first demonstrated for a prototype 
two-dimensional system. In Section HI we will define a "symmetric dual" of operators 
and states, to apply our formalism. In Section 5, the approach will be extended to the 
three-dimensional Bianchi I model and it will be shown that, with certain additional 
caveats, the same conclusions can be drawn. Finally, we will demonstrate that "lattice 
refinement" [HI |T6] can be motivated by alternative isotropic embeddings, while the 
difference from the standard quantisation of isotropic cosmology can be viewed simply 
clS db different choice of basic observables. 



2. Basics formalism of LQG/LQC 

The full LQG theory is uniquely derived from the requirements that the theory is 
diffeomorphism invariant and satisfies the SU(2) gauge freedom. In both, the full theory 
and in LQC, the fundamental variables are the holonomies of the SU(2) connection, A^, 
along a given edge, e, and the corresponding momentum, which turns out to be the flux 
of the triad through a two-surface S. The holonomies are given by 



hM) =Vexpj Y{s)A; (7(s)) t4s , (1) 

where V infers path ordering on the exponential, 7^ is the tangent vector along the edge 
e, and are the basis of the SU(2) Lie algebra. The corresponding momentum variable 
is 

E{SJ)= [ e^^.E'J.dx'^dx' , (2) 
Js 

where fi is an SU(2) valued test function and Ef is the densitised triad; k, - ■ ■ are 
SU(2) indices, whilst a,b,c, ■ ■ ■ are coordinate indices. 

Restricting to isotropic and homogeneous cosmologies allows us to consider only 
straight edges along integral curves of the basis vectors, Xf, that produce the network. 
Then the connection is given by a (dynamic) multiple of the basis one forms, o;^, namely 
= c(t)w*, whilst the triad is E^ = ^J^p{t)Xf , where is the determinant of 
the fiducial metricl§| and c stands for the connection component. This is simply a 
consequence of the symmetries imposed, since isotropy ensures that there can be no 
angular dependence of the connection or triads, whilst homogeneity ensures they are 
the same at every spatial point (but not necessarily on different time slices). 

With these symmetries, the holonomies become simply 



hi{A) = exp 



-«/ioO-i ~ 
c 



= cos i^^j + 2Ti sin l^^j , (3) 

§ The fiducial metric is a complication that arises only for an open universe and is used to define the 
volume to which spatial integrals are restricted to ensure they remain finite. Physical results should 
not depend on this volume |17j . 



Lattice Refining LQC from an Isotropic Embedding of Anisotropic Cosmology 4 



where r,, the basis of the SU(2) Lie algebra, are related to the Pauli matrices, cxj, by 
Tj = —iai/2 and /xq is the orientated length of the edge with respect to the fiducial 
metric. It is this /iq that is the ambiguity coming from the fact that we have imposed 
homogeneity on our system. In a full theory, the length of the edges would have some 
spectrum of values, in which case /iq would be the minimum such length, which has 
been shown to be non-zero [18]. Here we have treated constant, however for 

consistent quantisation [HI [9l |7], it has been shown that it must vary as /io ~ jj,"^/"^. 
In order to include this varying parameter, one must change variables from the triad 
component to the volume [5], however here we will briefly sketch the derivation for the 
constant /xq case, with more details given in Ref. [5]. 

By analogue with the full theory, the kinematic Hilbert space is extended via the 
Bohr compactification of the real line [18] . An orthonormal basis for this Hilbert space 
is {|/i)}, where 

{c\^^) = e^^ . (4) 

Note that c is a re-definition of the connection component, so that the Poisson bracket 
of the connection component and the (re-defined) triad component is independent of 
the volume of the fiducial cell, namely {c,p} = kj/S, with k, = SvrG and 7 the Barbero- 
Immirzi parameter, fixed through a black hole entropy calculation and turns out to be 
7 ^ 0.2375. 

The triad operator acts on these basis states as 

PW = -»^f |w = ^MW. (5) 

Clearly, also the volume operator, V = = has eigenstates in this basis, namely 



To calculate the eigenvalues of the inverse volume operator, the classical expression |18j . 

= , (7) 

is used. This is a classical identity independent of L, but will be quantised to different 
operators for different L. Thus, L plays the role of a quantisation ambiguity. 

It is now possible to quantise the classical Hamiltonian constraint, where the 
curvature is approximated via holonomies around a closed curve (which can not be 
shrunk to zero, due to the "area gap", given in terms of fiQ] it occurs by analogue with 
LQG), which gives 

Qrav = -i^tr Y: e'^' {h^^^h^^^h^^-'h^^-' \ht^^-\ V] sgn {p)) , (8) 

where, as in all quantum theories, there is an ambiguity in the factor ordering [9j; other 
quantum ambiguities have been fixed. Using Eq. (jS]), one can thus obtain a difference 
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equation [T7], 



V, 



+ 



V. 



/I— 5/10 



^m+4mo(</')-4 



A* -4/10 



(0) 



of step 4/xo, for the wave-function coefficients ^ (/x), defined by 



*/i(0) 



(9) 



(10) 



Note that, the matter Hamiltonian T-L^ is assumed to act diagonally on the basis states 
with eigenvalue Ti^^fi). 

By adapting the network to the symmetries of a Bianchi I model, one can similarly 
derive the quantum Hamiltonian for this anisotropic system |T3]. In this case, we 
have three triad components {pi,P2,P3), along the three directions of the anisotropic 
model. The gravitational sector of the kinematic Hilbert space consists of wave-functions 
^(Pi,P2,P3), satisfying 

^(Pi,P2,P3) = ^(bi|,b2|,b3|) • (11) 
Neglecting the details of the derivation of the Hamiltonian constraint, which are of no 
relevance for our present study, we draw the attention of the reader to the fact that in 
both, the isotropic and the Bianchi I case, the holonomies act on a state as shift vectors. 
Let us define three dimensionless variables, Aj, as 
sgn{pi)^/\pi 



^ ^ ^ v^.yv i^M j^j^ i = 1, 2, 3 ; (12) 

(47r|7|yA/|i)'/3 ^ ^ 

the quantum of area A/pj denotes the physical geometry, with A = An^y/S. The 

Hamiltonian constraint for the wave-function \E'(Ai, A2, i^), where z/ = 2A1A2A3, is a 
difference equation. It reads [1] 



TXG 



[z/ + 2)vV + 4^+(Ai,A2,z/; 

-(z/ + 2)757 ^o+(Ai,A2,z/; 4>) 
- (z/-2)v^^o(Ai,A2,i/; 4>) 
+ (z/ - 2)v/|z/-4|*j(Ai, A2, z/; 4 



(13) 



where \E'^4 are defined as [H] 

^^(Ai,A2,z/; 0) = 

^ ( Ai, 



V 



u±2 



+ ^ ( Ai,A2,z/;0 



z/±2 



■ Ai, • A2, l^, ' 



V 



+^ Ai 



V 



z/±2 



■ A2,z/;0 +^ A 



z/±2 
z/±2 



■ Ai, A2, v\ < 



^1, 



■ A2, z/; 



(14) 
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and 



z/±4 

V 



, z/±2 , 
Ai, A2, ± 4; : 



Ai, — — - ■ A2,z/±4;. 



z/±2 



z/±4 
z/±2 
z/±2 



Ai, A2, z/ ± 4; < 
• Ai,A2,z^±4;. 



Ai, A2,z/±4; 



+ ^ Ai 



z/±4 
z/±2 



A2,z^±4;. 



(15) 



respectively. Note that, plays the role of internal time, exactly as in the isotropic case. 
Thus, as one can easily check from Eq. (fT3|) . as far as the u dependence is concerned, 
the steps are uniform, since the argument of the wave-function involves z/ — 4, z/, z/ + 4, 
exactly as in the isotropic case. However, the dependence on Ai, A2 is quite complicated, 
technically, to deal with. 

Writing the step-size in terms of the triad components, implies 



K 



^iso (P) ^ ^iso P ± 



and ^ {pi,P2,P3) < 



^ [pi±K 



P2P3 ' 



^[PUP2±KJ\J^^ 



P2,P3 
P3 



(16) 



^ [Pl,P2,P3 ± K 



P3 
P1P2 



for some suitable constants K and K that are proportional to "vA/pi [H]. Note that, 
typically the K constants, appearing in the argument of the anisotropic wave-function, 
will be different. 

In both, the isotropic and anisotropic cases, "consistency" is derived entirely within 
the model, i.e., the symmetries of the model are fixed from the outset and "consistency" 
is looked for, under these assumptions. This, of course, does not ensure that the model 
remains "consistent" , if it derived from a less symmetric approximation. In particular, 
the isotropic case is consistently quantised for the step-size of the holonomies given in 
Eq. (|T6|) . only if one does not allow additional degree of freedom (and similarly for 
the Bianchi I model), however this does not ensure that an isotropic embedding of the 
Bianchi I model will have the same step-size as the isotropic case, given in Eq. (fT6|) . 
However, with both models fully quantised, it is possible to consider such an isotropic 
embedding of the anisotropic Bianchi I model and explicitly evaluate the type of step-size 
that is induced on the effective isotropic system. 

This will be done in some detail below, but one can immediately see that setting 
Pi = P2 = Pz = P, will produce the expected form of step-size, however this embedding 
procedure does not lead to a viable isotropic theory, instead one has to use a projection 
of the anisotropic degrees of freedom |13]. The reason for this is the following: for 
Pi = P2 = Ps = Pi we have Ai = A2 = A3 = {v/2Y^^ . Thus, the state which satisfies 
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Pl=P2=P3=PiS 




From Eq. (IT^ . we see that applying the Hamiltonian to such a state involves states of 
the form ^ ^Ai, A2, t'j , which do not satisfy Ai = A2 = A3 = (z//2)^/^. Thus, states for 
which Pi = P2 = Ps = P, do not form a super-selection of states (i.e., the Hamiltonian 
does not preserve this property of states). 

The purpose of this paper is to show that it is, in fact, possible to find a super- 
selection of states that are all (measured to be) isotropic. In particular, we will show that 
there are embeddings for which the expectation values of the variables pi, p2 and p^ are 
equal, and that are preserved under the action of the Hamiltonian. The Hamiltonian 
applied to such states does not produce the standard FLRW Hamiltonian of LQC, 
rather one finds a Hamiltonian constraint in which the discreteness length varies. Thus, 
one may view this as a form of lattice refinement on the isotropic sub-space. It is 
important to note however that the typical forms of lattice refining Hamiltonian used 
in the literature are not the same as those produced here; there is instead a simplified 
attempt to capture the missing underlying degrees of freedom (in this case the underlying 
anisotropic degrees of freedom) . 

3. The two-dimensional case 

To demonstrate the approach, let us first consider an anisotropic wave-function that 
depends only on two directions, labelled 1,2. The basic operators that are quantised 
are pi and hi, with i = 1,2. Note that pi = af, while hi denote the holonomies. 
In general, the wave-function is anisotropic, namely 



l^(Pl,P2)) 7^ |^(P2,Pl)) • 

To be clear, the action of the operators pi is 



(17) 



Pll* {Pl,P2)) = Pll^ {Pl,P2)) 
P2|^ iPl,P2)) = P2|^ iPl,P2)) 



(18) 



and 



Pll^ {P2,Pl)) 
P2|^ (P2,Pl)) 



P2|^ {P2,P1)) 

(P2,Pl)) , 



(19) 



and the action of the holonomies hi on a states |\E' {P11P2)) and |\E' {p2,Pi)) are, 

hl'^ iPl,P2)) = 1^ (Pi + 5i iPl,P2) ,P2)) 

and /i2|^ (Pi,P2)) = I* {pi,P2 + S2 (pi,P2))) , 



(20) 
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and 

hi\-^ {p2,Pi)) = 1^ iP2 + Si ipi,P2) 
and /lal* (P2,Pi)) = I* (P2,Pi + ^2 (Pi,P2))) , 

where 6i is potentially a function of both pi and p2- 



(21) 



3.1. Standard isotropic embedding 

Typically, one can go to the isotropic case by defining the operator Qsym, as 

4ym|^ iPl,P2)) = 1^ iPuPl)) = 1^ iP2,P2)) = |^ {p,p)) ■ (22) 

In other words, the operator Qsym sets pi = p2- It is easy to see that operating singly 
with either hi, or /i2, on (pi,P2)), before operating with Qsym, implies 6i = 0, e.g., 

Qsym hil"^ {pi,P2)) = Qisol^ iPl+ Sl,P2)) 

= |^(pi + + (5i)) 

= l*(P2,P2)) 

Pi + Si=p2 , (23) 

at least for singly defined wave-functions. Thus, the operation is not conserved under 
the action of a single holonomy, unless 6i = 0. This can be understood easily, since 
a general anisotropic state can be always made isotropic by setting pi = p2, but if we 
change either pi or p2, then the condition required to now make them equal will, of 
course, change. 

If we consider only isotropic applications of the holonomies, e.g., hih2, then the 
operation can be conserved under the action of Qsym, provided 6i = 62- To be explicit. 



sym 



-- <5sym|^ (Pl + Sl,P2 + 62)) 
= |^(pi +5i)) 
= |*(P2 + 52,P2 + ^2)) , 
Pl+Sl=P2 + 62 , (24) 

which preserves the condition pi = p2 iff 61 = 62, for singly defined wave-functions. 

So, we found that the value of p in either direction is the same and it changes in 
the same way, namely it is isotropic. Notice that, for consistency, this requires that we 
only consider isotropic applications of the holonomy operators [e.g., the Hamiltonian) , 
which is physically acceptable since if our system is to end up isotropic, then it must 
evolve isotropically. However, in the Bianchi Hamiltonian (Eq. (fT3l) ) the holonomies act 
via multiplication, rather than addition, of the variables thus. 



sym 



/ll/l2|^ (pi,P2)) 



5lPl 



= <5sym|^ {SiPi,62P2)) 

= |^(5ipi,(5ipi)) 

= I* {hP2,hP2)) , 

S2P2 , 



(25) 
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which again preserves the condition pi = p2 iff ^1 = ^2, (again for singly defined wave- 
functions). The full Bianchi Hamiltonian (Eq. f lT3|) ) contains linear combinations of 
operations of the form of Eq. ( l25l) that do not satisfy 61 = 62- Thus, the condition 
Pi = P2 is not preserved. In order for the standard "improved dynamics" of the isotropic 
model to be recovered, one must project out the anisotropic degrees of freedom jl4j . 



3.2. Alternative isotropic embedding 



There is also another approach one could follow to forming an isotropic system from the 
fully (two-dimensional) anisotropic state |\E' (^1,^2))- Let us define the symmetrisation 
operator by 



4ym|^ (Pl,P2)) = - [1^ {Pl,P2)) + |^ {P2,Pl)) 



(26) 



Thus, whilst (for example) pi Qsyml^ (P15P2)) 



does not appear to be isotropic, the 



physically expectation value of pi between two such states is. To be clear, 

1 
2 



Pi 



Qsyml^ {Pl,P2)) 



Pl\^ iPl,P2)) +P2\^ iP2,Pl)) 



(27) 



is not symmetric, however the inner product, 

( ^ (Pl , P2 ) Qsym |Pl I Qsym ^ (Pl , P2 ) ) 



1 
4 

(Pl +P2) 



Pi {Pl,P2) 1^ {P2,Pl)) +P2(^ {P2,Pl) 1^ {Pl,P2)) +P2 



[1 + (^(P1,P2)|^(P2,P1))] 



= (^(Pl,P2)Qsym|P2|Qsym*(Pl,P2)) , 

is clearly symmetric. Note that, we used {X\Y) = {Y\X) and the notation 



{^ipi,P2)Q: 



sym I 



Qsyml^ {P1,P2)) 



i2t 



(29) 



Thus, the measured scale factor along either direction is the same. 

We can now consider what the action of an holonomy is and, in particular, whether 
it is conserved under the symmetrisation operator. Consider, for example, the operation 
of hi. 

^ V (Pi + ^uP2)) + 1^ {P2,Pl + Si)) 



Q: 



sym 



hi\^{pi,P2)) 



(30) 



We can again check that this state is symmetric, in that the expectation values of the 
two Pi's are the same, namely 

(pi,P2) hiQsjm\Pl\Qsymhl'^ {Pl,P2)) 

P1+P2 + Si 



1 + {^{pi+6i,P2)\^{p2,Pl + Si)) 



ipi,P2) hiQsym\P2\Qsymhl^ iPl,P2)) ■ 



(31) 
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The above conclusion is however only valid, if the symmetrisation operator is always 
the last one to be taken on the state. This has some undesirable features, in particular, 

0(4ym|^)) 7^4ym(0|^)) , (32) 

namely that the symmetrisation operation does not commute with other operators, 
which is problematic, since the operation of a generic operator (O) on a symmetric 
state, does not, on its own, lead to a symmetric state. In the next section, we will 
show how that this problem is avoided provided we are interested only in "symmetric 
operators" . In order to define such operators, it is useful to define a process we refer 
to as the "symmetric dual" . 



4. The "symmetric dual" 



Here, we will define the "symmetric dual" of operators and states. Note that, the 
symmetrised version of an operator must be used even if pi = p2 (we refer the reader 
to the previous section), since one has to ensure that only symmetric products of the 
holonomies are taken. 



Define the symmetrised version of a state, |\l/)sym, as 



1^) 



sym 



^1 
^1 



*(Pl,P2)) + |^(P2,P1))] , 



(33) 
(34) 



where by |\I')^ we denote the symmetric "dual" of a state (pi,P2))- Similarly, define 
the symmetrised version of an operator, Osym, as 

1 



Osym — 



+ 0' 



where the operation of the "dual" operator is defined by requiring, 
0|^))^ = OS|^)S , 



and noting that (|\I^)^)' 



Csyml^') 



|\E'), by definition. We thus get 
1 



sym 



AB 
1 

A 



where 



1^) 



1 

B L 



6|^) + 6^1^) 



(35) 
(36) 

(37) 
(38) 

(39) 



Thus, we find that the operation of a symmetric operator on a symmetric state is (as 
expected) a symmetric state. This can be visualised on the two-dimensional (^1,^2) 
plane. The case of pi = p2 is the diagonal line, and if we restrict ourselves to states that 
lie on this line, then we must also restrict to operators that act diagonally, as one can 
see in Fig. [H However, if we consider a symmetric superposition of states, then in order 
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P2 




■ Pi = P2 



0|V) 



pi 



Figure 1. Visualisation of the symmetrisation of states. 



for the operation of a general operator O to give another symmetric state, one should 
require the operation to be symmetric. One easy way of doing this, is to define the dual 
operator, O^, and act with the symmetric combination Osym ~ {0 + 0^). Note that, 
this procedure is by no means unique, however such symmetric operators maintain the 
symmetric subspace made up of states |\l/)sym, i.e., if we restrict ourselves to considering 
only such operators, then we have a super-selection of states |^')sym- 

In Fig. [H we have taken the operator O to act on a state by shifting (^1,^2)- This 
is, of course, motivated by the fact that we want to consider holonomies, which we 
will assume act as shift operators, as in Eq. ( I20l) . Consider the restrictions placed on 
^i{Pi,P2) and 62 {pi,P2), by requiring that 

hih2 = 6 = 6^ , (40) 

i.e., that the action of a holonomy along each direction, is a symmetric operation in the 
sense defined above. As we will see, does not mean that the state is "shifted" diagonally. 
We have 



hlh2\'^) = /il/i2|^ (Pl,P2)) = 1^ (Pl + Si (P1,P2) ,P2 + S2 (Pl,P2))) 



and 



/il/i2|^)^ = hih2\'^ {P2.P1)) = 1^ {P2 + Si {p2,Pl) ,Pl + S2 {P2,P1))) 

Thus, to satisfy 

r "1 S 

one should require that 

1^ iP2 + ^2 iPuP2) ,Pl + Sl (Pl,P2))) = 

1^ {P2 + Si {p2,Pl) ,Pl + S2 {P2,P1))) , 



(41) 



(42) 
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which imphes 

^2 (Pl,P2) = ^1 iP2,Pl) ■ (43) 

For illustration, let us consider some simple examples. For 6i = 62 = const., the pair of 
holonomies acts diagonally, with a constant step. A more interesting example is 

1 1 



Si 



5o 



which is valid only when pi = P2, namely on the diagonal line. This is the analogue 
of the "new" quantisation (or "improved" quantisation) approach. Motivated by the 
three-dimensional anisotropic results of Ref . [11] , one can take 



Pi 



^1 {Pi,P2) = \ — 

P2 



S2 iPl,P2] 




Clearly there are many other cases one may consider. 

The question then is, what would be the measured step-size between states 
separated by this symmetric operator. We have 

sym(^ (Pl,P2) (Pl,P2))sym = sym(^ (Pl,P2) |P2|^ (Pl,P2))sym 

P1+P2 



P1+P2 

A^ 
P1+P2 



l + {^{pi,P2) \^{P2,P1)) 
1 + (^1^)^ 



(44) 



where in the last line, we used the fact that the symmetrised state should be normalised, 
and we have chosen 



(^I^)S = (p^^p2) I* (P2,Pl)) = A^-1. 

Defining, 

\lp (Pl,P2))sym = hih2\'^ (pi,P2))sym , 

and taking the inner product, we find 

sym(V' iPuP2) |Pl | ^ (Pl , P2) )sym =sym (Pl,P2) |P2 | (Pl , P2) )sym 
P1+P2 + Si (pi,P2) + 5l {P2,Pl 



(45) 

(46) 



A^ 



1 + {^{pi,p2) ^(P2,P1) 



(47) 



where 
and 



iPl,P2) = {Pl + ^1 (Pl,P2) ,^2 + ^2 (Pl,P2)) 



(P2,Pl) = {P2 + ^1 {P2,PI) ,Pl + ^2 {P2,Pl)) ■ 

We can then evaluate what the measured step-size, Agym, would be: 

Asym = sym(V^ (^1,^2) (Pl,P2))sym -sym (^1,^2) \Pi\^ {Pl,P2))sym 

= {pi+p2){X -l) + {6i{pi,p2)+S2{pi,P2))X , (4^ 
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where X is defined as 
Thus, we have 

1 + (^l^)s ' ^ ^ 

and we used \x) = \^ {pi,P2)) to emphasise the possibihty of a difference from 
|^(Pi)P2))- The normalisation of the symmetrised state gives the denominator (see, 
Eq. fH5|) ). but we get X = 1 iff the norm of every state (or every state in the super- 
selection, if one is possible) is the same. In particular, if we were to restrict ourselves to 
Pi = P2 = p, this would imply that pi = p2 and hence that \x) = \x)^ and |\^) = |\^)^. 
This gives X = 1, which means A = 5i(pi,P2) + {P2,Pi) = 25 (p). However, for 
general ("symmetric dual") states, X can be greater, or less, than unity. Thus, the 
step-size would include a different dependence on pi and p2] it can be seen as a "lattice 
refinement" analogue. 

In Section [5l we will extend to three directions and explicitly use the Bianchi I 
model of Ref. [H] to evaluate the step-sizes for the two different symmetrisation schemes. 
Indeed, the above can be used directly for the case that two of the directions are set 
equal to each other. For example, if the restriction Pi = P2 ^ Ps is made, then we 
are restricted to a sub-class of Bianchi I models, which are essentially two-dimensional. 
Such cylindrically symmetric models would be useful for investigating black holes, where 
the expectation value of p is different along the radial and angular directions. 

Intuitively, what we have is that the use of superpositions to symmetrise a state, 
means that (for symmetric operators, such as pairs of holonomies with the correct type 
of step) we can find symmetric states which however have different step-sizes than in 
the standard case. Essentially, the standard case is just a special case of the above 
procedure, when the two states used in the superposition are the same. 

To be clear, we have chosen to normalise the states, such that 

{^{puP2)\^{p2,Pi))=A'-l. (51) 

However, the restriction 61 (^1,^2) = ^2 {P2,Pi), is not sufficient to simultaneously specify 
the normalisation of 

[Pl + 61 {pi,P2) ,P2 + S2 (Pl,P2)] 1^ [P2 + Si (p2,Pl) ,Pl + ^2 (p2,Pl)]) • 

For Pl = P2, 

Asym = 5i {pi,P2) + §2 iPl,P2) 

= 26{p) , (52) 

where the last equality follows from the fact that 61 (^1,^2) = ^2(^2, Pi) = ^{p)- For 
any other choice of (^1,^2), such that pi ^ P2, this is not the case. Whilst we cannot 
rule out the possibility that the (symmetric) state evolves towards a state made up of a 
superposition of two even more anisotropic states, this seems unlikely. 
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If we take that, at least for large pi + P2, we have (pi, J92)) — |\1' (p, p)), then 
we expect X \ 1 {i.e., the inner product of the "symmetric dual" states tends to 1 
from above; it can be understood by thinking of the states as vectors in a vector space). 
Thus, 

Asym > 5l (j5l,J92) + (^2 (j5l,J92) , (53) 

whilst tending to equality as pi and p2 become large and tend to each other. 

For the case of an isotropic system, made up of ever more anisotropic states, the 
converse is true, namely, 

Asym < ^1 (Pl,P2) + '^2 (Pl,P2) , (54) 

whilst not necessarily tending to equality as pi and p2 become large. Note that here we 
also require Eq. ( I43l) to hold, which is what ensures that the measured step-size, Agym, 
is indeed symmetric under interchange of pi,P2- 

In conclusion, in general the step-size can be either larger or smaller (and has a 
different dependence on (^1,^2)) than in the pi = p2 case. For the case of a state 
made up of a superposition of ever more isotropic states, the measured step-size tends 
to the standard step-size as pi, p2 become large. This is precisely what is modelled 
by "lattice refinement "[iJ], albeit in an heuristic way, except that we might expect the 
"lattice refinement" to approach the special case of the "new" quantisation refinement 
rate, for large scales. 

5. The three-dimensional case 

In this section, we demonstrate a concrete example of the above procedure for the case 
of an anisotropic Bianchi I model, recently consistently derived in Ref. [H]. In this case 
the states are defined on the orthonormal basis |Ai, A2, A3), where Aj are essentially the 
scale factors along the three directions of the Bianchi model. 

In the standard approach, isotropic states are taken to be those in which the 
dependence on Ai,A2 and A3 is "integrated out", or more precisely, defining the volume 
of the state as z/ = 2A1A2A3 to eliminate one of the directions (A3), the map 

\^ (Ai, A2, z/)) ^ I 5^ VI/ (Ai, X2,iy)) = \^il (z/)) , (55) 

produces isotropic states. This projection produces exactly isotropic system with the 
equations given in terms of the volume (called the "z/ quantisation"). Here we will show 
that it is also possible to have a super-selection of states that are isotropic and which 
evolve according to an alternative quantisation procedure which may be modelled by 
lattice refinement. We are going to work with the three scale factors A,, to demonstrate 
that there is an ambiguity in exactly what the measured volume of such isotropic states 

II Here and throughout we use "lattice refinement" to refer to any quantisation in which the lattice 
defined with respect to /i is dynamically varying. In particular, then, "new" quantisation is simply a 
special case of such "lattice refinement" . 
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would be. Firstly, we want to find a symmetric superposition of states that is isotropic, 
in the sense defined in Section I3.2[ 
Consider the state 

l^r (Ai, A2, A3)) = ^ [1^ (Ai, A2, A3))+|^ (A3, Ai, A2)) + |^ (A2, A3, Ai))] , (56) 

with the additional restriction 

(^(Ai,A2,A3)|^(A3,Ai,A2)) 
= (^(Ai,A2,A3)|^(A2,A3,Ai)) 
= (^(A3,Ai,A2)|^(A2,A3,Ai)) 

= / , (57) 

on the anisotropic states being used. This additional requirement is not present in the 
two-dimensional case, since there is only one cross-correlator between two anisotropic 
states, which is trivially equal to itself. In the three-dimensional system however, 
discussed here, we have three cross-correlators, which need to be equal if the state 
is to be isotropic. This is a restriction that would be present in any dimension higher 
than 2. 

The expectation values of the scale factors along each direction of such a state are 
/\ \ Ai + A2 + A3 

(Ai) = • (58) 

Note that we have chosen A to normalise the state. It is worth mentioning that, by 
choosing a different state, the conditions required for the expectation values to match 
change. In particular, for the state 

\X (Ai, A2, A3)) = J [{^ (A3, Ai, A2) 1^ (A2, A3, Ai)) 1^ (Ai, A2, A3)) 
+ (Ai, A2, A3) 1^ (A3, Ai, A2))|^ (A2, A3, Ai)) 
+ (Ai, A2, A3) 1^ (A2, A3, Ai))|^ (A3, Ai, A2))] , (59) 

to be isotropic, we require 

|(vl/(Ai,A2,A3)|vI>(A3,Ai,A2))P 

= |(vl/(Ai,A2,A3)|vi/(A2,A3,Ai))|2 

= |(vI/(A3,Ai,A2)|vl/(A2,A3,Ai))|2. (60) 

With such a state, one also finds that (Aj) = (Ai + A2 + A3) /3, which explicitly 
demonstrates the non-uniqueness of the choice of isotropic embedding. 

For both these states, the measured scale factor is equal in each direction and is 
given by the average of the scale factors of the underlying, anisotropic states. However, 
the measured volume of such a state, found using the volume operator, 

(Ai, A2, A3)) = 2AiA2A3|^ (Ai, A2, A3)) , (61) 

is just u = 2A1A2A3, which is not the cube of the measured scale factor. Thus, whilst it is 
the eigenvalue of the anisotropically defined volume operator, it is not necessarily what 
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we would measure as the volume. Essentially, this is because whilst the scale factors A, 
are measured to be equal in each direction, they are not eigenvalues of the state, i.e., 
Ai|^ (Ai, A2, A3)) 7^ Ai|\^ (Ai, A2, A3)). However, both the average and the product (i.e., 
the volume) of the scale factors are eigenvalues. It is this ambiguity, that leads to the 
possibility of deviations from the standard isotropic case. In particular, for our choice 
of isotropic states, we have essentially defined the projection, 

|^(Ai,A2,A3))^^((A),z/((A))) , (62) 

where the state depends only on the average, (A), of the three scale factors, and their 
product {i.e., the volume) can be viewed as a complicated function of this average. 

The Hamiltonian found in Ref. [11] is linear and can be written as a symmetric 
combination of three operators that act on each Ai, A2 and A3 equally. Thus the 
Hamiltonian operator satisfies our definition of a "symmetric operator", Eq. ( !35l) and 
hence we can directly apply the above states to it. Alternatively one can directly evaluate 
the Hamiltonian acting on a symmetric state and show show that it contains only states 
that are themselves symmetric. 

As it is explained in Ref. [1^, the Hamiltonian acts on uniform steps along the 
u direction, just as in the "new quantisation" of isotropic cosmologies. Just as in the 
standard case, one finds that the Hamiltonian, acting on symmetric states, produces a 
super-selection of points along the u direction, that are equally spaced, however from 
Eq. (l62l) it is clear that this no longer implies a uniform spacing in either (A) (the 
measured scale factor) or (A)^ (the measured volume). To deduce the step-size in the 
measured scale factor (which would give the step-size in the holonomies that should be 
used in the case of isotropic LQC), one would have to invert the function z/((A)). The 
use of "lattice refinement" in isotropic LQC is a heuristic first approximation to this, 
and whilst the power law like refinements typically used in the literature |8j are rather 
crude, they have the advantage of being analytically tractable. 

There are many ways that such isotropic states can be produced from superpositions 
of anisotropic states. Certainly, those given in Eqs. (156!) and (1591) are not the 
most general ones, however they serve to demonstrate the possibility that whilst the 
underlying, anisotropic system requires specific, well defined holonomy step-sizes, the 
isotropic version contains an ambiguity in the precise value of the step-size that one can 
take. Here, we have focused only on the Bianchi I model, for which a full understanding 
of the quantisation exists, however one expects similar effects to occur when starting 
from other anisotropic and inhomogeneous models. In particular, the motivation for 
lattice refinement in LQC comes from the fact that one expects that the additional 
degrees of freedom present in the underlying, full theory, will not exactly produce the 
symmetry reduced models of cosmology. Here we have shown that the set of states 
that are measured to be isotropic form a super-selection of states within the Bianchi 1 
model, which have a non-constant measured volume step. In this sense we have shown 
that the anisotropic degrees of freedom of the Bianchi 1 model, that are missing in the 
symmetric FLRW case, can produce a non-standard evolution of the FLRW volume. 
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which is exactly what is modelled by lattice refinement. It is important to note however 
that the Hamiltonian governing the evolution of these isotropic states is not the standard 
FLRW Hamiltonian that is usually considered in lattice refinement models. Thus whilst 
it may be hoped that lattice refinement models capture some of the consequences of a 
varying discreteness size, the use of the standard Hamiltonian makes this a rather crude 
approach. 

Finally, it is important to note that in Eq. fl62|) . we chose to view the measured scale 
factor of the isotropic state as the primary variable, however one could have chosen to 
view the eigenvalue of the volume operator as fundamental, in which case the isotropic 
volume of the state would be u and the isotropic Hamiltonian would indeed produce a 
super-selection with uniform steps in volume. In this case however the scale factor of the 
state ((A)) would be related to the volume through a complicated function (essentially 
the inverse of i/((A))). This would seem to be largely a matter of taste, however it is 
possible that future work may choose one view point over the other. 

6. Conclusions 

Historically, the quantisation of FLRW cosmologies within LQC, considered the triad 
element to be the fundamental variable that is to be quantised |18]. However, it was 
realised that this leads to an instability within the theory at large scales [H], which is 
particularly problematic for infiation [7]. It was later shown that these problems are 
cured [7] for the "new" quantisation approach in which the physical volume is considered 
as the fundamental object to be quantised [5]. 

In terms of the basic triad component p, the discretisation of the Hamiltonian 
occurs with step-size given by p"^, where ^ = for the original quantisation approach 
and A= —1/2 for the "new" quantisation. These two possibilities are the two limits in 
which, either only the labels of the edge spins, or the number of vertices dynamically 
change in the underlying lattice (respectively for the "old" and "new" quantisation 
approaches). Whilst the behaviour of the full Hamiltonian constraint on a general 
lattice has yet to be derived from LQC, it is expected that the constraint on A should 
be < ^ < —1/2, since the Hamiltonian acts by a combination of both processes. A full 
derivation of this from an underlying LQC point of view, is not yet possible, and here 
we do not restrict the value of A and refer to the general case as "lattice refinement" El- 
There has been much work (and significant debate) on whether "lattice refinement" for 
a general A is theoretically consistent O El [9] and strictly working within the symmetry 
reduced FLRW models, it appears that only for ^ = — 1/2 ("new" quantisation) is this 
the case [6]. 

% Although the onset of inflation within LQC may remain problematic, even with the use of the "new" 
quantisation [19]; a question which has been recently addressed in Ref. [20j. 

+ The term "lattice refinement" is also sometimes used to include only the cases with —1/2 < A < 
case, which explicitly excluded the "new" quantisation approach. Here, as previously, we use the term 
"lattice refinement" in its most general sense, to refer all quantisations in which the step-size is given 
by p-^, which includes both the "new" and "old" quantisations as special cases. 
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However, here we have shown that, given a (fully consistent) quantum anisotropic 
model, it is possible to find isotropic states for which the discrete step-size of the 
isotropically embedded Hamiltonian constraint is not (necessarily) that of "new" 
quantisation. This demonstrates the possibility of "lattice refinement" being due to 
the degrees of freedom that are absent in the isotropic model. It is of course possible to 
embeco the "new" quantisation approach within the anisotropic models and it may be 
that this is the more natural procedure, however without an a priori knowledge of which 
embedding process to follow, the conservative approach would be to consider both. 

We have shown that the difference between the two procedures is essentially due 
to what one considers to be more fundamental, the volume of the underlying states 
(z/) or the measured volume of the symmetric state ((A)^), which are not necessarily 
equal. Choosing u leads to the "new" quantised Hamiltonian of isotropic cosmology, 
whilst choosing (A)^ results in some kind of different "lattice refinement". This "lattice 
refinement" is significantly more complicated that the single power law behaviour, p-^, 
usually considered, however one may hope to capture some of its effects using the simple 
model. 

It is important to note that here we have started from a Bianchi I model and looked 
for isotropic embeddings, however similar ideas should apply more general models. 
We are not suggesting that the Bianchi I model is the underlying geometry of our 
universe, however it is sufficient to demonstrate the ambiguities associated with finding 
an effective isotropic cosmology within a model that has more degrees of freedom. This 
can be the thought of as the source of the "lattice refinement" ambiguity of LQC, 
however more work is needed to find out if the simple power law model is a sufficiently 
accurate approximation of this ambiguity. More work is also necessary to see how this 
new formulation of "lattice refinement" , in particular the new relationship between the 
effective scale factor, (A), and the volume, u, effects the consistency of the isotropic 
theory. 

Finally, one may be concerned with the fact that simply deciding which variables 
are the more fundamental one, can have important implications for the large scale 
physics of the theory, however this is not an unusual occurrence for quantum systems. 
In all quantum systems, we start from a classical action and follow specific procedures 
to produce the quantum version of the theory. Precisely which classically equivalent 
action to choose is a source of quantum ambiguities. From a theoretical point of view, 
one may hope to have the choice picked out from a more complete understanding of 
the underlying theory, whilst on a practical level, it is often more useful to parametrise 
our ignorance and fix these parameters from experiments and observations. In the 
latter case one simply assumes that the approach is an effective theory, approximating 
some, as yet unknown, full underlying theory, which is sufficiently accurate for the 
current data. In the case of LQC (and indeed classical cosmology) we are constantly 
working within the framework of an effective theory of universal dynamics, in which (it 

* Here the term "embedding" is used in a loose sense. The FRLW model is in fact a "projection" of a 
system with more degrees of freedom. 



Lattice Refining LQC from an Isotropic Embedding of Anisotropic Cosmology 19 



appears) the symmetry reduction is typically a valid approximation. In this context, 
"lattice refinement" is simply an effective parametrisation of the fact that cosmological 
symmetries are not exact symmetries of the physical universe. Whilst measuring any 
consequences of this effective parametrisation seems unlikely, it is possible that the 
largest scales of the Cosmic Microwave Background (which correspond to the earliest 
scales to leave the horizon) might retain a signature of this LQC regime, from which 
ambiguities such as "lattice refinement" may be fixed. 
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